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Abstract. We introduce a class of Kac-like kinetic equations on the real line, with general 
random collisional rules, which include as particular cases models for wealth redistribution in 
an agent-based market [6_, or models for granular gases with a background heat bath 
Conditions on these collisional rules which guarantee both the existence and uniqueness of 
equilibrium profiles and their main properties are found. We show that the characterization 
of these stationary solutions is of independent interest, since the same profiles are shown to be 
solutions of different evolution problems, both in the econophysics context 6 , and in the kinetic 
theory of rarefied gases 1141 1291 . 



1. Introduction 

In this paper, we are concerned with the study of the time evolution and the asymptotic behavior 
of the spatially homogeneous kinetic equation 

j dt^t + IH = Q + (Ht,Ht) 

1 M(0) = 



(1) 



which caricatures a Boltzmann-like equation in one spatial dimension. The solution fi t = A*t( - ) is 
a time-dependent probability measure on R, describing, in its most common physical applications, 
the distribution of particle velocity in a homogeneous gas, which is initially distributed according 
to the probability measure /2o- The gain operator Q + models velocity changes due to binary 
particle collisions. Our fundamental assumption is that Q + is a generalized Wild convolution. 
More precisely, for all bounded and continuous test functions g £ Cf,(R), we characterize the 
probability measure by 

(2) J g(v)Q+(n,fi)(dv) =E[JJg(y 1 A 1 + v 2 A 2 + A Q )fi(dv 1 )fi(dv 2 ) 

where (Aq, Ai, A2) is a random vector of M 3 defined on a probability space (f2, T, P) and E denotes 
the expectation with respect to P. 

The interaction rule generated by the law described in ([2]) simulates an interaction in which, 
in addition to the standard binary collision, the post-interaction velocities are randomly modified 
by the presence of an external background. As we shall see, this modification induces an evolu- 
tion process for the probability measure which stabilizes in time towards a steady profile heavily 
dependent of this random collision part. The physical relevance of this generalized collision rule 
is mainly related to the dissipative Boltzmann equation. Indeed, in a dissipative binary collision 
process, a particular choice of this random contribution is shown to produce the same steady state 
of the classical Boltzmann equation with standard dissipative binary collisions, in presence of a 
thermal bath [14] . 

A second main example of application of equation @ is linked to the field of econophysics [B] . 
In this case, the generalized collision refers to a market based on binary trades between agents, in 
which part of the traded money is taken away by an external third subject, which redistributes it 
according to a certain economical random rule. 
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For A = and for suitable choices of (A 1; A 2 ), the one-dimensional kinetic equation ([T]) 
reduces to well-known simplified models for a spatially homogeneous gas, in which particles move 
only in one spatial direction. The basic assumption is that particles change their velocities only 
because of binary collisions. When two particles collide, then their velocities change from v and 
w, respectively, to 

(3) v = p\v + qiw w' — P2V + q2W 

where {p\,qi) and {q2,P2) are two identically distributed random vectors (not necessarily inde- 
pendent) with the same law of (Ai,A 2 ). 

The first model of the type ([T])&([2]) has been introduced by Kac [22], with the collisional 
parameters Pi — sin9 and % = cos#, i = 1,2, for a random angle 9, uniformly distributed on 
[0, 27r). The dynamics describes a gas in which the colliding molecules exchange a random fraction 
of their kinetic energies. This idea has been extended in [25] to gases with inelastically colliding 
molecules, which loose a random part of their energy in each interaction. The inelastic Kac 
equation corresponds to ([T])&([2]) with pi = |sin6>| p sin# and qi = | cos#| p cos#, with p > being 
the parameter of inelasticity. Recently, more general versions of ([T])&([2]) have been considered: 
their applications range from gases under the influence of a background heat bath to models for 
the redistribution of wealth in simple market economies 18, 123] . In most of the above mentioned 
cases, A\ and A2 are positive random variables such that K[Af + A|] = 1 (conservation of energy) 
[221 [15], or +A2] = 1 (conservation of momentum) [15] 123], 

In the classical Boltzmann equation [I2J [13] relaxation to Maxwellian equilibrium (Gaussian 
density) is shown to be a universal behavior of the solution. Contrary, the corresponding equilibria 
of model (TTJ), to which the solution is shown to relax, depend heavily on the precise form of the 
microscopic interactions ([3]). Furthermore, they are not always explicitly known analytically. 

In the case of models of wealth distribution in the society, the comparison of these steady 
states with realistic data is up to now the only means to evaluate — a posteriori — the quality 
of a proposed model. For instance, it is commonly accepted that the wealth distribution should 
approach a stationary (or, in general, a self-similar) profile for large times, and that the latter 
should exhibit a Pareto tail [16l [17]. The asymptotic behavior of the solutions of ([T]), when 
Aq — 0, has been extensively treated in [2] [3[ [9], and it is by now fully understood, in particular 
if one aims to describe a few analytically accessible properties (e.g. moments and smoothness). 

The general situation in which Aq ^ 0, while relevant in various applications which will be 
dealt with in this paper, has never been touched before. This case corresponds to assume that in 
a binary interaction the particle velocities change from v and w, respectively, to 

(4) v' = piv + q\w + 771 w' = p 2 v + q 2 w + r/2 

where (pi , q\ , 771 ) and (92 , Pi , 772) are two identically distributed random vectors with the same law 
of (Ai, A2, Aq). We will now describe the specific examples we are dealing with. 

Kinetic models of a simple market economy with redistribution. In [6] Boltzmann-type kinetic 
models for wealth redistribution in a simple market economy have been introduced and discussed. 
The authors focused their attention to models which include taxes to each trading process. As- 
suming that a percentage of the total wealth involved in the trade is not returned to agents, the 
goal in [B] was to understand the role of redistribution, there produced by a linear transport-drift 
type operator. 

Here we assume that the economic trades between agents are described by an interaction of type 
(j4|). In particular, for a given positive constant < e < 1, the post-interaction wealths including 
redistribution are given by 

(5) Ai = (1 — e)A 1 A 2 = (1 - e)A 2 A Q = eA 

where M[A% + A2] = 1 and E[Ao] = mo = J vp,o(dv). Note that within this assumption, the total 
mean wealth is left unchanged. The mixing parameters (Ai, A2) can be chosen among the variety 
of models present in the pertinent literature, see e.g. [23 . The classical model introduced in [15) 
corresponds to the choice 

Ax =X + f,{l-X), I 2 = rj(l-A) 
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where rj is a random variable denned on [0, 1] (symmetrically distributed around 1/2) and A £ [0, 1] 
is a parameter (the so called saving propensity), while the pure gambling [5J corresponds to fix 
Ai=A 2 = fj. 

An interesting variant of the previous model is obtained by setting 

(6) Ai = (l- eA)Ai A 2 = (l- eA)i 2 A = eAi 

where {A\, A 2 , Aq) and A are stochastically independent, and P{A = 1} = 1 — P{A = 0} = 8. 
The presence of A in (|6|) simulates a market in which taxation does not act on the totality of 
trades, but it occurs only with a probability 6. 

As we shall see in Section [321 one can nx the values of (A±,A 2 ,Ao) in such a way that the 
steady state of the model (p} produces the same steady states as the model considered in [6] . 

Inelastic Kac models with background. A second interesting application of binary interactions of 
type (|3]) is related to the study of a dissipative gas in a thermal bath [TTJ [T3]. In one space- 
dimension, a dissipative Kac-like model has been introduced and discussed in [35J- As already 
mentioned, this model corresponds to the choice 

(7) A x = |sin((9)| p sin(<9), A 2 = | cos{8)\ p cos(0) 

where 9 is uniformly distributed on [0, 2ir). As shown in 4, 25], in consequence of the dissipation, 
a solution to the Kac equation corresponding to an initial value with finite second moment con- 
verges in time toward the probability mass located in zero. In addition to the physical dissipative 
interaction ([7]). let us now assume that particles velocities are subject to random fluctuations rji, 
induced by an external background, whose distribution is the same of Aq, while Aq and (Ai,A 2 ) 
are stochastically independent. In addition let us assume that Aq ^ 0, but E[Aq] = 0. 

As extensively discussed in |11) , and directly verifiable on the single binary collision, the presence 
of this random fluctuation of zero mean is such that the post-collision energy is bigger than the 
corresponding one induced by the dissipative collision without fluctuations, i.e. when Aq = 0. 
Indeed, since A and (Ai, A 2 ) are stochastically independent and EL4 ] = 0, 

E {{v'f + {w'f) = E(A 2 + A 2 )(v 2 + w 2 ) + AEiA^vw + 2E(A 2 ) 

= E(A lV + A 2 w) 2 + E(A lW + A 2 v) 2 + 2E(^) 

with E(Aq) > 0. The main consequence of this fact is that one can exhibit examples in which the 
initial value has finite second moment and at the same time the corresponding solution does not 
converge in time toward a degenerate distribution. The same phenomenon is shown to happen if 
one adds to the dissipative Boltzmann equation a thermal bath [14] . 

This allows to establish a direct link between the steady states of the present dissipative colli- 
sional models with random fluctuations and the steady states of the dissipative Boltzmann equation 
in presence of diffusion [14], as well as in presence of friction and/or drift [29]. Indeed, the steady 
states of the various problems on the dissipative Boltzmann equation quoted above, are steady 
states of the Boltzmann problem ([T]), corresponding to suitable choices of the random variables 
(Ai, A 2 , Aq). We will detail the correspondences between these problems in Section [3] 

2. Main results 

We start by writing the Boltzmann equation ^ in Fourier variables. By setting tp(t, £) = 
/ e l ^ v Ht{dv) , and using Bobylev's identity [7], one obtains that </>(£,£) obeys to the equation 




d t <p(t, o + 4>(t, e) = Q+ (<M*, -),<Kt, •)) (0 (t > 0, £ £ R) 



where 

(9) q+(<kv),<Kv))(0 -^n^A^^A^y^ }. 

The initial condition </>o(£) = / e^ v p,Q(dv) can be seen as the characteristic function of a prescribed 
real random variable Xq, i.e. </>o(£) = E[e l ^ x °]. 
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As in the case of the Kac equation, it is easy to see that (|5J| admits a unique solution (f> which 
can be written as a Wild series [3T] 

(10) 0(t,o = E e ~ t ( 1 - e_t ) n 9«(o. 

n>0 

where <7o(£) = <A)(£) an d, for n > 1, 

71—1 

(11) 9«(0 = -E ( 3 + fe^n-l-.)(0- 

Hence, if /it is the unique solution of (flj with initial condition /2q, then its Fourier- Stieltjes 
transform is given by (|10p . 



2.1. Steady states. The stationary equation associated to (jSJ is 

(12) 0oo(O=<3 + (0oo,0oo)(O K6R). 

It can be proven that, under suitable hypotheses, a solution to (|12[1 exist. To show that steady 
states exist it is enough to recast the problem as a problem of fixed point equation for distributions. 
In terms of probability distributions, ([T2")l reads 

(13) Q+(M, M )=/i, 

where, given any probability distribution //, by ([2]), the probability distribution Q + (p,,fi) is the 
law of the random variable 

A + YxAi + Y 2 A 2 , 

Yx and Y2 having law \i and Y\, Y2 and (Aq, A\, A2) being stochastically independent. 
In what follows, let us set 

My ■= probability measure on £>(R) : J \xp fi(dx) < +00 1 , 

and, for every m in K and 7 > 1, 

My^m := j(i£ A4 7 : / xfx(dx) = m 

Finally, when E[Ai + A 2 ] ^ 1, let us define 

E[4>] 



' 1 - EL4 X + A 2 ] ' 

The convex function q : [0, 00) — > [0, 00] defined by 

(14) g( 7 ):=E[|A| 7 + |A 2 H, 

where 0° := 0, will play a very important role in what follows. 

First of all, let us collect some known results on the existence of solutions of equation (|13p . 



Proposition 2.1 ([27], [28]). Assume that there is 7 in (0,2] such that E[|Ao| 7 ] < +00 and 

q(i) < 1. 

(a) If < 7 < 1, i/ien i/iere is a unique solution /i^ 0/ (| 1 3[) m My. In addition, if 7 = 1, 
t/iis solution belongs to M\.m,; 

(b) 7/ 1 < 7 < 2 and EL4i + A2] 7^ 1, then there is a unique solution fioo of (| 13() in My and 
this solution belongs to My,rh> 

(c) // 1 < 7 < 2, E[Ai + j4a] = 1 and E[j4o] = 0, then, for every mo € R, there is a unique 
solution /Xoo 0/ (Ti"5)) m .M 7 , mo . 

Let us notice that in case (a) (cfr. Lemma [5.2[1 , it is possible to describe fi^ in terms of a 
suitable series of random variables. 

While it is easy to check when is a degenerate distribution, necessary and sufficient con- 
ditions for boundedness of moments up to a certain order are more difficult to obtain. A partial 
answer to this problem is given in the next proposition. 
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Proposition 2.2. Let the same hypotheses of Provosition \KT\ be in force. 

(i) In case (a) or (b) of Proposition 0O1 fi^ is a degenerate distribution if and only if m(l — 
(Ai + A%)) — Aq almost surely (a.s.) for some real number to; in case (c) of Proposition 
\2.1\ /Uqq is a degenerate distribution if and only if mo (I — [A\ + A2)) — Aq a.s.; 

(ii) If q((3) < 1 and E[|t4o|^] < +00 for some f3 > 2, then q(s) < 1 for every 7 < s < f3 and 
J \x\ fj,oo(dx) < +00; 

(iii) Let Aq, Ai and A 2 be positive random variables with P{Aq ^ 0} > 0. //, for some 
f3 > max{l,7} ; J \x\@ ^^(dx) < +00 and J x^i OD (dx) > 0, then /ioo{[0, +00)} = 1 and 
qifi) < 1. 

2.2. Trend to equilibrium. We recall that the Kantorovich-Wasserstein distance of order 7 > 
between two probability measures and v is defined by 

(15) Uu,u):= inf (E\X' - Y'\'<) 1 ' max ^) 

' (X',Y') 

where the infimum is taken over all pairs {X 1 , Y') of real random variables whose marginal prob- 
ability distributions are fj, and v, respectively. 

If iy n ) n is a sequence of probability measures belonging to M. 1 and G A / f 7 , then l^(u n , Voo) — » 
as n — > +00 if and only if v n converges weakly to and 

J \x\"'v n {dx) -> J \x\~tvooidx). 

See, e.g., [26] . Recall that v n converges weakly to means that J g(x)u n (dx) —¥ J g[x)v 00 [dx) 
for every g in Cb(M). 

We are now ready to state our main results concerning the long time behavior of the solutions. 

Proposition 2.3. Let 7 G (0, 1). Assume that E[|Xo| 7 + |A)| 7 ] < +00 and (7(7) < 1. Let fi^ be 
the unique solution in A4~ ( to (fl"3f . Then, for every t > 

l-y(th, Poo) < ^(Mo,Moo)e~* (1_9(7)) . 

In what follows, whenever E|Xo| < +00, set mo = E[Xq]. 

Proposition 2.4. Assume that K[\Xq\ + \Aq\] < +00 and that q(l) < 1. Let fi^ be the unique 
probability measure in Mi which satisfies (|13[) . Then for every t > 

and J v^j, 00 (dv) = fh. Moreover, if mo = fh, then J v^it(dv) — fh for all t > 0. 

Proposition 2.5. Assume that, for some 7 G (1,2], E[|X | 7 + \A \~<] < +00, 5(7) < 1, E[A 1 + 
A%] =/= 1 and to = to. Let fj.^ be the unique solution in A4 7 to (|13[) . Then, for every t > 0, 

« 7 (M t ,Moc) < 2 1 / 7 L r ( Mo , Moo )e- t ( 1 -^»/ 7 
and f R v[ioo(dv) = fh. Moreover, J vfit(dv) — fh for all t > 0. 

Proposition 2.6. Assume that, for some 7 G (1,2], E[|Xo| 7 + |^4o| 7 ] < +00, 5(7) < 1 and that 
E[j4i + A2] = 1 and M[Aq] = 0. Let fj,^ be the unique solution in .M 7;mo of (|13[) . Then 

Z 7 (Mt,Moo) < 2 1 / 7 ; 7 ( Al0 ,Moo)e- t(1 - 9(7))/7 
for every t > 0. Moreover, J vfit(dv) = too for all t > 0. 

3. Examples 

3.1. Kinetic models of a simple market economy with redistribution. The first applica- 
tion of the results of Section [2] deals with the kinetic model for wealth with redistribution, briefly 
described in the Introduction. In this leading example, the random variables A\,Ai, Aq are given 
by ([6]). As already noticed, these assumptions correspond to a kinetic model for wealth distribu- 
tion in which part of the wealth put into the binary trade is taken away by a third subject, which 



6 



FEDERICO BASSETTI, LUCIA LADELLI, AND GIUSEPPE TOSCANI 



at the same time restitutes to agents a certain amount of wealth. This is done in such a way that 
the mean total amount of wealth into the system is left unchanged. Assuming ((5]), one has 

E[Ai + A 2 ] = 1 - e5 < 1. 

Hence, since E[Ao] = / vp,o(dv) = mo < +oo, one can invoke Propositions Ell] (a) and !2.4l to prove 
both the existence and uniqueness in A4i of a steady state and the (exponential) convergence to 
this steady state of any solution with finite initial moment of order one. 

One of the interesting effect of the redistribution is that the steady state can have finite moments 
of higher order than those of the steady states of the corresponding model without redistribution. 
This can be easily verified by comparing the steady states corresponding to e > 0, say (j!oa , with 
the steady states without redistribution, we will denote by obtained by setting S = e = 0. 

Thanks to Theorem 5.3 in [2D], it is known that, given /3 > 1, J v@ fj!^ \dv) < +oo if and only 
if q(/3) := E[Af + A$] < 1. On the other hand, if e > and E[A$] < +oo, by Proposition EH 
(ii)-(iii), / v p n£ e] '(d«) < +oo if and only if q(fi) < 1. Since q(/3) = [1 + S[(l - ef ~ l]]g(/3) 
and [1 + S[(l — c)P — 1]] < 1, one can easily give examples in which J (dv) = +oo while 

fvPfi&^idv) < +oo. 

The previous discussion does not solve another interesting problem connected with wealth 
taxation and redistribution: the existence of an optimal amount of taxation. If one assumes that, 
given a certain (conserved) amount of money, the optimal redistribution refers to a steady state in 
which all people in the market ends up with almost the same amount of money, this problem can 
be solved by looking for the steady state with minimal variance. We leave this point to a further 
research. 

3.2. Connections with other form of redistribution. We show here that the law of (A\, A2,Aq) 
can be fixed in such a way that the steady states of the redistribution model proposed in [6] fit 
into our framework. 

Let us start by briefly outlining the model introduced in [6j. In Fourier variables this model 
reads 

(16) o + <P(t, o = m, o + Kim, o 

with 

(17) Q £ (^4>):=E[ct>{At^(AU)}, 
and 

(18) R £ x (m) ■■= -ext-^m + ie(x + 1)"*>W(0 (x > -1), 

In p7]) (AIjA^) are positive random variables such that E[A\ + AJj] = 1 — e, and -g?(f>(Q,Q) = 
i J vp,o(dv) — im . Note that in (|16[) the interaction operator consists in a dissipative collision 
operator, given by Q e (<j),(j)), and a redistribution (differential) operator i?^.(</>)(£). It is worth 
recalling that, if <f> is the Fourier- Stieltjes transform of a (regular) density /, then Rt(4>) is the 
Fourier-Stieltjes transform of 

(XV- (x + l)m Q ) /(«) . 
The possible steady states of (fTB|) must satisfy 

(19) 0(O = 4(0,0)(O+^W)(O- 

Existence of a global solution t) to ([16]) has been proved in [6] provided that J vflo(dv) — mo. 
Anything was proven about the existence (and eventually uniqueness) of a steady state. This 
problem can be solved in a surprisingly easy way by establishing a connection between the steady 
states of the model [DJ and special cases of our model. 
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First of all let us fix \ = —1 in ([Taj) . In this case the redistribution operator simplifies, and 
equation (flU]) reduces to 

(20) m =&(<M)(£)+e^(£). 

Resorting to the analogous computation in Bobylev, Cercignani and Gamba [HI [H] , equality (|20l) 
can be equivalently rewritten as 

(21) <}>{£) = / Q^,m^ e )du. 



Jo 

It is immediate to see that equation (12T1) can be rephrased as 

(22) cfr(O=E[cl>(U-*A* 1 Z)cl>(U- e A*0}, 

where U and (Al,A 2 ) are stochastically independent and U is uniformly distributed on [0,1]. 
Hence, the steady state (j2"0j) coincides with the steady state (JT2J) corresponding to (Ao, Ai, A 2 ) := 
(0, U~ e A\, U~ e A 2 ). Since in this case q(l) = EL4i + A 2 ] = 1, in order to apply Proposition 12.11 
(c) it is necessary that A\ and A 2 satisfy 

= ^fUlT + [Alf] < l, 

for some 1 < 7 < min{2,l/e}. If E[(Aj;) 7 + (Aj) 7 ] < 1 this inequality holds true for every 
e < (1 - E[(Aj) 7 + (A^)' 1 })/^. It should be noticed that in this case, since A = 0, and Ai and A 2 
are positive with EL4i + A 2 ] = 1 one can resort also to Theorem 2(a) of |20) . 

Let us now consider the case in which \ = 0, so that (|18[) corresponds to a pure transport 
operator, which produces a uniform redistribution. In this case, equation (|19[) becomes 

(23) 0(O = &(0.0)(O+femo#(O. 
or, what is the same, 

(24) m = i-t mo A ^,m)- 

Let us observe that if Aq is an exponential random variable of mean emo, that is with density 
ho(v) — exp{-u/(erao)}/(emo) (v > 0), then 

E[e* Ao ] = / e* v ko(v)dv = -. 

Jo 1 - lemot; 

Under the additional assumption that Aq and {A\,A 2 ) are stochastically independent, ([24]) can 
be equivalently written as 

^)=E[e^<j>{A\^{Al()]. 
Hence, it is enough to choose A\ = A\,A 2 — A 2 and Aq as above to identify the steady state (|23l) 
with the steady state lfl~2|) . Note that, since in this case EL4i + A 2 ] = 1 — e < 1, the assumptions 
of Proposition [23] (a) are trivially satisfied for 7=1. 

Last, let us examine the physically relevant case in which x > — 1 and x / 0. For any given 
e G (0, 1] set S := e\ > ~ e - With this choice, (fT§)) becomes 

4W>, M) = <K0 + K-j^m - i(S + e)mo^(e)- 
Multiplying both sides for e _, ^ mo ^ i ^ _1 we get 



_ d_ 

which, integrating over [0,£], gives 



ft 

(25) 6e- iimo ^^(f)^)= Qei^^i^e'"" 10 ^^- 1 ^. 
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By the change of variable £u — r, we can write the previous equation in the equivalent form 

»i 

iu 



Jo 

which can be rephrased as 

0(0 = E[</)(C/ 5 A*O0(^^Oe l(1 ' t/ ' ) ^ mo5 ] ; 

where (A^A^) and U are stochastically independent and U is uniformly distributed on [0,1]. 
Hence p^|) is equivalent to 

0(0 - EMA^MAzOe^} 

for 

(26) A 1 = U 5 A* 1 A 2 = U S A* 2 A = (1 - U s )^-^m . 

o 

Since 5 > — e, it is immediate to reckon that E[A\ + A 2 ] = (1 — e)/(l +5) < 1. Hence Proposition 
12.11 (a) applies. In particular, since E[v4o]/(l — E[Ai + A2]) = mo, the solution described in 
Proposition 12.11 satisfies f v[ioa{dv) = niQ. 

3.3. Inelastic Kac models with background and connection with dissipative models 
with diffusion. A further application of the results of Section [21 announced in the Introduction, 
results from the choice 

(27) A 1 = |sin((9)| p sin09), A 2 = | cos(0)| p cos(6>) 

with 9 uniformly distributed on [0, 2tt). This assumption leads, when Aq = 0, to the inelastic 
Kac model [25] , which describes the cooling of a one-dimensional spatially homogeneous Maxwell- 
like gas. In fact, if A = and J \v\ 2 ^ 1+p ^ fj,o(dv) < +00, fi t is shown to converge weakly to the 
probability mass concentrated in (cfr. [4]). As we shall see, the addition of a random fluctuation, 
described by the random variable Aq ^ 0, is responsible for the formation of non-trivial steady 
states. If A\ and A 2 are given by (|27|) . then 

E[Ai + A 2 \ = 0, 

and, whenever 7 > 2/(1 + p), 

9(7) = E[|Ai| 7 + \A 2 p] = A- f\\ S in(0)| (p+1)7 + I cos{e)\^ +1 ^]de < 1. 

27r Jo 

Hence, if for some positive e, with e + 2/(1 + p) < 2, E|A | 2/(1+p)+£ < +00 (with E[A ] = m := 
/ v/j, (dv) if p < 1) and J \v\ 2 ^ 1+p ^ +€ ft (dv) < +00, then Propositions 12 .311231 apply. In particular, 
if E[A ] = and P{A ^ 0} > 0, the steady state is a non-degenerate probability distribution 
with finite moments of all orders. 

As a special case let us choose A = A , a — with ^o,a and A ,b exponentially distributed 
with density v n> exp{— v/a}/a and v h-> exp{— v/b}/b (v > 0), and assume that Ao ta , Ao t b, Ai, A 2 
are stochastically independent. Since 

1 1 



Efe^ ] = R[e* A °-']E[e 



1 (l-io£)(l + i6f) l-i{a-b)t + ab£?' 
if a := (mo + \J m\ + 4c 2 )/2 and b := (— rno + \/ m o + 4cr 2 )/2 the stationary equation (|12[) becomes 

E^xO^O] - 0(0(1 - im Q i + <r 2 e) 
which can be equivalently written, after setting Qp(<fi, 0) := E[cf>(Ait;)<p(A 2 £)], as 

(28) 0(o = q+(0, 0)(o - ° 2 em + im o e0(o- 

Equation (|28|) describes the steady states of the inelastic Kac equation in presence of a thermal 
bath and a transport term. Indeed, if is the Fourier- Stieltjes transform of a density /, then 
— c 2 £ 2 0(£) + imo£0(£) is the Fourier-Stieltjes transform of 

<9 2 9 
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In particular, the analysis of Section[2J allows to prove existence of a steady state for the dissipative 
Kac equation with diffusion. The problem of the solvability of equations of type 

(29) Q(/,/) + <r 2 A/ = 0, 

in terms of nonnegative integrable densities / £ L^R 3 ), and where Q is the Boltzmann collision 
operator, is a well-known problem in kinetic theory of rarefied gases. When Q is the dissipative 
collision operator for Maxwellian molecules, existence of non trivial weak solutions has been proved 
by Cercignani, Illner and Stoica [14] . 

Also, as clearly discussed by Villani in [29], apart from collisions, other physically relevant 
problem in kinetic theory of granular gases lead to the addition of various terms which either model 
external physical forces, or arise from particular situations. One of these situations is described 
by equation ([25]) . A second one is obtained by subtracting a drift term to the Boltzmann collision 
operator. This leads to the problem of finding steady states of the equation 

(30) Q(fJ)-a 2 V-(vf)=0. 

Let us remark that in one dimension of the velocity space, equation (|30[) is a particular case of 
equation (|19j) with ^ = — 1, which has been solved in the previous Sub-section. 

4. Probabilistic representation of the solutions 

The core of the proofs of our results is a suitable probabilistic representation of the solution 
fit ■ The idea to represent the solutions of the Kac equation in a probabilistic way dates back, at 
least, to the work of McKean [23], but it has been fully formalized and employed in the derivation 
of analytic results for the Kac equation only in the last decade, starting from [TU] and [2Tj . 

Our approach here follows the same steps used in [2] and [3] and it is based on the concept of 
random recursive binary trees. It is worth recalling that a binary tree is a (planar and rooted) 
tree where each node is either a leaf (that is, it has no successor) or it has 2 successors. We define 
the size of the binary tree r, in symbol |r|, by the number of internal nodes. Hence, any binary 
tree with 2k + 1 nodes has size k and possesses k + 1 leaves. Any binary tree can be seen as a 
subset of 

U = {0}U[U fc >i{l,2}*]. 

As usual is the root and if a = (<Ji, ... ,ak) {&i £ {1:2}) is a node of a binary tree then 
the length of a is \a\ := k. Moreover (a, (Jk+i) ■— (cij • • • 1 0>, Cfe+i) and for every 1 < i < k, 
a\i :— (<7i, . . . , <Ji) and er|0 = 0. 

We now describe a tree evolution process which gives rise to the so called "random binary 
recursive tree". The evolution process starts with To, an empty tree, with just an external node 
(the root). The first step in the growth process is to replace this external node by an internal one 
with 2 successors that are leave. In this way one obtains T\. Then with probability 1/2 (i.e. one 
over the number of leaves) one of these 2 leaves is selected and again replaced by an internal node 
with 2 successors. One continues along the same rules. At every time k, a binary tree with k 
internal nodes is obtained. For more details on binary recursive trees see, for instance, [19] , 

In the rest of the paper, given a binary tree r, we shall denote by £(r) the set of the leaves of 
r and by X(t) the set of the internal nodes of r. 

The Wild series expansion (|10l) - (fTT| can be translated in a probabilistic representation of the 
solutions as sums of random variables indexed by binary recursive random trees. On a sufficiently 
large probability space (f2, J 7 , P) let the following be given: 

• a family (X v ) ve u of independent random variables with common probability distribution 
Mo; 

• a family (Ao(v), Ai(v), ^2(^)) ugU of independent positive random vectors with the same 
distribution of (A , A\, A 2 ); 

• a sequence of binary recursive random trees (T ra ) ng N; 

• a stochastic process (vt)t>o with values in No such that P{v t = k} = e~'(l — e~ t ) k for 
every integer k > 0. 



in 
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Write A(v) = (A (v), A\{v), A 2 {v)) and assume further that 

(A(v)) V £ V , (T„) n >i, {X v ) veV and (vt)t>o 

are stochastically independent. 

For each node v = (v±, . . . , Vk) in U set 

\v\-l 



i=0 



and w(%) = 1. Define 
and, for any n > 1, 



W := X $ and T := 



W n := J2 ™( v ) x v> r„:= E ro H^o(«), W*:=W„+r w . 
we£(T„) vei(T n ) 

Proposition 4.1. Equation |8]) /ias a unique solution <j>, which coincides with the characteristic 
function of W* t , i. e. 

oo 

4>(t,0 = E[e 1 ^] = e~*(l - e-*) n E[e iew »] (t > 0, £ e R). 

n=Q 

Proof. We need some preliminary results on recursive binary trees. A very important issue is 
that any binary tree has a recursive structure. More precisely we can use the following recursive 
definition of binary trees: a binary tree r is either just an external node or an internal node with 2 
subtrees, r^ 1 ' , r^ 2 ' , that are again binary trees. For every k > let denote the set of all binary 
trees with size k. By Proposition 3.1 in [2], we know that if (Tk)k>o is a sequence of random 
binary recursive trees, then for every k > 1, j = 0, . . . ,k — 1 and every r in Tfe, 



(31) 

and for k > 1 
(32) 



,{ T (i) =T (i) )T ( 2 ) =T ( 2 )| |T ( 1 ) |=i } 



= P{Tj = r^jPin.j., = t W}1{\t™\ = j} 



= t( 2 )1MI t < 1 )| = 



for every j = 0, . . . , k — 1. Now observe that, in order to prove the proposition we need only to 
prove that q n (0 = Eje^^j, for every n > 0. This is clearly true for n — 0. For n > 1, write 



\v\-l 



w: = Aom+j2 A M{[ E Ii A vlMi)xi j) ] + [ E II 41(^)4^)]} 

J = 1 i>e£(T«>) 1=0 veX(Ti 3) ) i=0 

where = A(0»), A<r° = X {j<v) and, by convention, if C(T£) = the ter ms between 

square brackets is equal to = Xj. Since (A^(v), Xy )„etr, j — 1> 2, are independent, with the 
same distribution of (A(v), X„)„eu> using pip and the induction hypothesis one proves that 

2 

(33) E [e^\A(<D), |T«|, \T^\] = J[ g |T u )| (^ i (0))e^W. 

3=1 

At this stage the conclusion follows easily by using (|3"2"j) : indeed: 

2 1 n-l 

E[e^*] =E[j]^ )| (eA J -(0))e^W] = - £ E[«&(£A 1 )g„_ J -_ 1 fcA 2 )e**>] = 5n (£). 



3=1 



3=0 



□ 
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5. Proofs of Section [2] 

Proposition 12. 1[ for 7 = 2, is proved in [37], while for general 7 G (0,2) it can be seen as a 
special case of a (more general) result contained in [28]. The proof of Proposition 12 . II given in [28] 
is based on some contraction properties of the Wasserstein metrics. Here we provide a proof for 
7 G (0, 2] based on a martingale method inspired by [27] . In this way we obtain some additional 
information on the solution, used to prove Proposition 12.21 

In the following we need to consider a sequence (T*)„>o of (deterministic) binary trees. Any 
such tree can be seen as a subset of U: starting from Tq = 0, for each n denote by T* the binary 
tree obtained from T*_ ± replacing each leaf by an internal node with 2 successors. Recall that 
C{T*) (1(T*), respectively) denotes the set of the leaves (the internal nodes, respectively) of T*. 

Define K(/i) := Q + (/i, /i) and K n (S m ) as the n-iterate of the transfomration K applied to the 
mass probability concentrated on the real value m. Finally set 

M*:=m J2 Yl m ( v )Mv)- 

In the rest of the paper U 1 will stand for L 7 (f2, J 7 , P). 

Lemma 5.1. Let 5(7) < 1 for some 7 in (1, 2] and E|Ao| 7 < +00. Assume either 

(i) E(Ai +A 2 )^1 and m = E(A )/(1 - E(Ai + A 2 )) = m or 

(ii) E(Ax + A 2 ) = I, E(A ) = and m arbitrary. 

Then 

(a) ie l (S m ) is the law of M*; 

(b) (M*)„>o is a martingale with respect to {Q* L ) n >\, with Q* = a(A(v) : v G T*_ 1 ), such that 
E(M*) = m for every n; 

(c) sup n E|Af*| 7 < +00, hence (M*)„>o converges a.s. and in L 1 to a random variable 
such that E(M^) = m and E|Af^| 7 < +00; 

(d) the law /ioo of is a solution of ()13[) in 

(e) // (ii) holds true then 

M^^mZoc + Y^ w(v)A {v) 

n>0v£C(T») 

where Z^ is the almost sure limit of^ v£C ^ T ,^ vo{v) for n — > +00. 
Proof, (a) is immediate for n = 1. In fact 

K\S m ) = V(A 1 (9)m + A 2 (<D)m + A o (0)) = £>( ^ w(v)m + ro ( v )A)(v)J> 

v£C(T*) v£X(T*) 

where, for every random variable X, T>(X) denotes the law of X . Now, by induction, we obtain 
K n (S m )=K(K n - 1 (5 m )) 

v<EC(T*_ 1 ) v£T(T*_ 1 ) 

At this stage, denote by T* 1 (T* 2 , respectively) the left (right, respectively) binary subtree of T*. 
For every v in U and i = 1,2 set m l (v) = w{{i,v))/Ai(%) if A ( (0) ^ and w\v) = if A t (%) = 0. 
It is plain to check that 

Y vj l {v)m+ Y rol (w)A ((i,w)) 
(34) ve£(T*i) vei( T n) 

i = 1,2 are independent random variables with the same law of M*_ x . 
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Hence, 

K"(5„ 1 )=P(A 1 (0)( ]T w\v)m+ Y, ™\v)A ((l,v))) 

( 35 ) +a 2 (0)( X! ^»™ + E ^»A)((2,«)))+a>c 

As for (b) is concerned, clearly M* is integrable and Q* t measurable. Moreover, 
M* = M*_, + Yl ™(v)[(Mv) + Mv) - l)m + A (v)} 

and hence 

e[m;\q*_ 1 ]=m*_ 1 

in both cases (i) and (ii). Furthermore, E[M*] = m for every n. Since M* is a martingale and 
1 < 7 < 2, we can apply the Topchii-Vatutin inequality, see e.g. [T], to get 

n 

e[|m*h < e[|m *h + 2^e[|m; -m;_x| 7 ] 

n 

= m 7 + 2^E[| ^ m{v)[{A 1 {v)+A 2 {v)-l)m + A (v)} 1 . 

Now, since E[(Ai(w)+A 2 (w)-l)m + A (?;)|g*_ 1 ] = E[(Ai(t>) + A 2 (u) - l)m + A (u)] =0, by the 
Bhaar-Esseen inequality (see [30]) we obtain 

E[|M*H <m 7 + X]TE[ ^ K«)n 

3=1 vECiT^J 

where K = 4E[|(Ai + A 2 — l)m + A) I 7 ] < +°o by assumption. Now it is easy to see that, for every 
k > 0, 

E[ Y \Mv)r]=q(j) k 
veC{T*) 

with q(j) < 1. Hence, sup n E|M*| 7 < +oo and, from the elementary martingale theory, it follows 
that (M*) n >o converges a.s. and in L 1 to a random variable such that E[Af^] = m and 
E[|M^| 7 ] < +oo. The proof of (c) is completed. In order to prove (d) set 0„(£) = E[exp(i£M*)]. 
By (|34l) . it is clear that ([55)1 is equivalent to 

MO = E[0„_ 1 (A 1 e)0„-i(^ 2 )e 4 « Ao ]. 

From (c) we know that 4> n (0 converges to </>oo(£) = E[exp(i£M^)] as n — > +oo. Hence, by 
dominated convergence theorem, we get 

MO = EIMAOM^b)^ ] 

and the proof of (d) is completed. Arguing as in the proof of (c) it is easy to see that under (ii), the 
terms J2 v ec(T') zv{v)m and J2 v &x(t*) ^i^Aoiv), which form M*, are both uniformly integrable 
martingales. Hence (e) follows easily. 

□ 

Lemma 5.2. Let 5(7) < 1 for some 7 in (0, 1] and E|Ao| 7 < +00. Then, for every m, 

(a) lC l (8 m ) is the law of M*; 

(b) J2vec(T*) w { v ) converges to in L~> ; 

(c) r* = X) 1 , e x(T*) w ( v )Ao(v) is a Cauchy sequence in L 1 and hence it converges in L 1 to 
the random variable 

T *oo = Y E ^)A,{v)- 

n>0v£C(T') 



KINETIC MODELS WITH RANDOMLY PERTURBED BINARY COLLISIONS 



13 



(d) M* converges to in L 1 and the law fi^ ofT^ is a solution of (|13[) in M.~; 



Proof. The proof of (a) is the same of the proof of (a) in Proposition 15.11 Furthermore, since 
7 6(0,1], 



E 



E 



< E 



E m«) 

vec(T*) 



IT 



9(7) n - 



Using the fact that 5(7) < 1, (b) follows. In order to prove (c) observe that, if n > m, 

n-l 



|r;-r^n = E 



I E E w ( y )Mv) 

j=mv£C(T') 



1 



<E[|AoH E 9(7)^0 

j=m-\-l 



for n,m —¥ +00, i.e. (r* )„ is a Cauchy sequence in L 1 . As for assertion (d) is concerned, combining 
(b) and (c) we obtain that (M*)„ converges in L 1 to T^. Finally, arguing as in the proof of (d) 

□ 



of Proposition 15 . 1 1 we obtain that T*^ is a solution of (|13j) in M 1 . 

Proof of Proposition [Ql The existence of a solution fi^ of (fl"3|) in A^ 7 , as required in (a), (b) 
and (c), is a consequence of Lemma l5.1l and Lemma l5.2l Let us prove the uniqueness. Let fj,± and 
fi2 be two solutions of (fT3")) in M.~ r Let (Yv) v and (Y« 2 ^)ti two sequences of independent random 
variables such that, for every ueD, Yj (Yj 2 \ respectively) has law /ii (/12, respectively) and, 
in addition, 

(Y}%, (Y v W) v , (A(v)) v 

are stochastically independent. Then, following the same lines of (a) in Lemma 15.11 and Lemma 
5.21 it is easy to see that 

E ™(v)Y v W+ E ™(v)A (v) 
«e£(r*) «ei(r*) 
has law ^t, (i = 1,2). As a consequence, if 7 6 (0, 1], then 

/ 7 (/i 1 ,M2)<E[| E n7(«)(K«-rJ 2 ))| 7 

<E[|r (1) -r (2) HE[ E M«)H 
= E[|r (1) -r (2) Hg(7r 

and 5(7)™ — > for n — > +00. Hence ^1 = /12 and this proves (a). As far as (b) is concerned notice 
that if 7 e (1, 2] and EL4i + A 2 ] ^ 1, it follows that 

E[Y«] = E[yj 2 )] = m 

and applying the Bhaar-Esseen inequality 



( Mi , M2 )<e[| e ^x^ (1) -n (2) )| 

veC(T*) 

<2E[|y (1) -r (2) HE[ E M y ) 



= 2E[|y (1) -y ^n g ( 7 r 

and hence [i\ = ^2 again. The case (c) follows in an analogous way since we need to consider only 
fi± and fi2 in A4 7 . m (i.e. we fix the mean). □ 

Proof of Provosition \2.SX The proof of (i) is straightforward. 

The proofs of (ii) and (iii) are inspired by the proof of Theorem 5.3 in |20j . Let us first prove (ii). 
Note that, since q is a convex function, for every A in [0, 1], g(A7 + (l — A)/?) < \q{-f) + (1 — \)q(/3), 



I 
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and hence q(s) < 1 for every 7 < s < (3. In addition E|^4 | s < +°° since E|ylo|' 9 < +00. Now fix 
s < (3, with l<fc<s<fc+l, k integer. Then, for Xi > 

(36) ( = (X->) ' < E^I + E^iaW^^ 3 )^ 



for suitable constants Cj 1 j 2 j 3 and ji are integers such that ji < k and j\ + j'2 + j/3 = k + 1. Using 
(f36| it is easy to see that 

(37) E[|y XJ 4i + F 2 A 2 + A | s ] < g(s)E|y| s + Cl E[\Y k ]i + c 2 

if Y, Y\, Yi are independent random variables with the same law v and (Y, Y±, Y2) is independent of 
(Ai, A 2 , ^4o)- The constants ci and c 2 may depend on (3 but not on za Obviously (|37[) is equivalent 
to 



(38) / |a;| s (Ki/)(da;) < g(s) / \x\ s v(dx) + c x \x\ k v(dx) 



C2 



Let us consider first the case in which 7 € (1, 2]. Choose either m = m if E[/li + A 2 ] ^ 1 or 
to = too = E[Xo] if E[yli + A2] = 1 and EL4o] = 0, and take v = S m . From Lemma \5. II we know 
that K"(5 m converges weakly to /i^ and that 



sup / |x| 7 (K™£ m )(dz;) < +00. 

n J 

Let us now prove that if for k > 1 and k < s < k + 1 , one has 

(39) sup / |a;| fe (K™<5 m )(da;) < +cx) and q(s) < 1 

n J 

then 

(40) sup J \x\ s (K n S m )(dx) < +00 and J \x\ s ^(dx) < +00. 

In fact, applying iteratively (|38l) starting from ^ = 5 m , since J |x| s <5 m ((ia;) = |to| s < +00, we 
obtain 

n — 1 

/ |x| s (K",5 m )(dx) < |TO| s g ( s )" + C^g(^ 

for a suitable constant C. Hence, since q(s) < 1, one gets sup n J |x| s (K™(5 m )(rfa;) < +00. Further- 
more define g M (x) = \x\ s I { \ x \< M} + M s I {lx \ >M} then 

\x\ s fiacldx) = / liminf gM^^idx) < liminf / g M (x)iJ loc (dx) 

J M— >+oo M->+oo J 

< liminf lim / g M (a;)(K n «5 m )(dx) < lim / |x| s (K™(5 m )(da;) < +00. 



Now, since 7 > 1 and /3 > 2, then (|39|) is true for k = 1 and s = 2. As a consequence we obtain 
([40)) for s = 2. Let us iterate this procedure for k < k with k < /3 < k + 1. The last step starts 
from the validity of (|39[) for k — k and s = (3 whcih implies J ^^/^^(da;) < +00. 
If < 7 < 1 then 

E\A + AiYx + A 2 Y 2 \ <E|Y|g(l)+E|A |. 
Since q(l) < 1 and E\A \ < +00 we get 

/n— 1 
|z|(K"£ m )(^) < |TO|g(l)"+C^g(l)^ 

and hence, thanks to Lemma [5 



sup / |x|(K n <5 m )(c?a;) < +00 and / \x\fioo(dx) < +00 



At this stage (|39f is proved for k = 1 and s = k + 1 = 2 and we can go on as in the previous case. 
This proves (ii). 
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Let us prove (iii). If 7 < 1, from Lemma [5.21 (c) one obtains that //^{[O, +00)} = 1. Now 
assume that 7 <E (1, 2]. Since P{A > 0} = 1 and P{A Q ^ 0} > 0, then E[A ] ^ and only case 
(b) of Proposition 12.11 has to be considered. By assumption rfi = J x^, 00 {dx) > and hence from 
Lemma [5TT1 (c)-(d), since M* is positive a.s. for every n > 1, again ^{[O, +00)}. 

As a consequence, using the fact that j3 > 1 we can write 

X? = v (A + AxX x + A 2 X 2 f > Aq + A{X{ + A%xg 

{= D denotes the identity in distribution) if X, X\,X 2 are independent random variables with law 
^oo and (X,Xi,X 2 ) and (Aq, A\, A%) are stochastically independent. Then 

E[Xf] > q(j3)E{X?] +E[Ao] > ?G9)E[*f] 
since we are assuming that P{Ao ^ 0} > 0. Hence q{(3) < 1. 

□ 

Let us state a useful result which is proved, with slightly different notation, in Lemma 2 of [3] 
(see also Proposition 4.1 in [2]). 

Lemma 5.3. Let 7 > such that 9(7) = E[|Ai| 7 + IA2I 7 ] < +00. Then, for every n > 0, 



ve£(T„) 



Proof. Given the sequence (T„)„>i of random binary recursive trees, one can define a sequence 
{V n )n>i of U-valued random variables such that 

T n+1 =T n U{(K ) l),(V n ,2)} 

for every n > 0, where Vo = and V n G C(T n ). The random variable corresponds to the 
random vertex chosen to generate T„ + i from T„. Hence, by construction, 

1 



P{V n = v\T u ...,T n } = I{veC(T n )} 



n + 1 



for every n > 1. Since T = and ro(0) = 1, E[X^s£(t ) l n7 ( u )l 7 ] — 1 an d bence (|4"Tj) is true for 
n = 0. For n. > 1, 



[ J2 m«)I 7 ]=e[ e K«)r[(iAi(«)r+|A2^)i 7 -i)i{K-i=v}+i 

ue£(T n ) ' «e£(T„_i) 

-e[ £ w , )P ]( 1 + «feL± 



E 



ue£(T n _i) 

from the independence assumptions and since (Ai(v), A 2 {v)) has the same law of (Ai,A 2 ) for 
every v. Hence, by induction, 



n 

<[ E i-(«)r]=n( 1 



E| 

' veC{T n ) 3=1 

and (|4Tj) is proved. 



g( 7 )-l\ r(«( 7 )+n) 



r(n + i)r(«( 7 )) 



□ 

In the following denote by the law of . 

Lemma 5.4. Assume that, for some 7 in (0,2], E[|Xo| 7 + |A)| 7 ] < +o° a71 ^ 9(7) < 1- 
(a) IfO < 7 < 1, tfien 

^(C)Moo) < c n (7); 7 (^ ,Moo) 

/or every n > 0, where /loo is unique solution of (j L3|) in A^ 7 . Furthermore, if 7 = 1 
and E(Xo) = m i/iera E(W*) = J vC,*(dv) = fh for every n > 1. 
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(b) If 1< 7 < 2, E(A 1 +A 2 )^1 and E(X ) = m, then 
(42) i?(C/*oo)<2c n (7)^(W),/ioo) 

where fj,^ is the unique solution of (113|) m A / ( 7 . Furthermore E(W*) = J v£*(dv) = m 
/or every n > 0. 

(c) J/ 1 < 7 < 2, E(A ) = 0, E(Ai + A 2 ) = 1, E(X ) = too (mo arbitrary), and ^ is i/ie 
unique solution of (fT3| m A4 7t m , then (|42"T) ZioWs and E(W 7 "*) = J v(*(dv) — too for 
every n > 0. 

Proof. The existence and uniqueness of in the three cases (a), (b) and (c) is guaranteed by 
Proposition 12.11 On a sufficiently large probability space (fi, J 7 , P) consider a sequence (Y„)„eu, 
such that 

• ( J 4(u)) t , eU ,(T n )„>o, and F t ,) 1 , 6U are independent; 

• (X V ,Y V ) are independent and identically distributed for v in U, and each (X V ,Y V ) is an 
optimal transport plan for / 7 (/2o 5 Moo)) i-e. the law of X v is /io, the law of Y v is /ioo and 

Let us set U* — J2 v <ec(t ) Y v w(v) + T n . We now show that, for every n, the law of U* is /t^. In 
point of fact 



CG£(T„_i) 

exp {W ^ o7(u)y„ + r n _i + ro(u)( Ai(^i + A 2 (u)y e2 + jog*) )) } 

— d Y- 

= e[ ^ I{^ n _i = «}, 



se£(T„_i) 

where V n is defined as in the proof of Lemma 15.31 Hence, by induction, U* has the same law of 
lg, that is /Zoo. Now denote by Q the er-field er(A(u) : v S U, (T„)„>i) and observe that 

^ (1 ' 7) (CMoo) < E|F< - C/:p =E[E[| £ w{v)(X v - Y v 

uS£(T„) 

<fc 7 E[ W{v)\~<e\\X v -Y v 

v<£C(T„) 

The last inequality is immediate for 7 < 1 with k~ t = 1 while, if 1 < 7 < 2, it follows with k~ t — 2 
from Bhaar-Esseen inequality, since E[X„] = E[3^] which implies E X)ue£(T ) w i v )(X v ~ Y v ) 
0. Hence, using (jUJ), 

^(i, 7)(c * )/ioo) < \m(v) |t]/--(i.t) (/i0j ^ = fc 7 c„( 7 )^ ax(1 ' 7) (Mo^oo). 

ue£(T„) 

In order to conclude the proof, let us study E[W*] when 7 belongs to [1,2]. Observe that 

E[W^]=E[W:_ 1 ]+e[ ™{v)[A (v) + A 1 (v)X vl + A 2 (v)X v2 -X v ]I{Vn = v} 

uer(T n _i) 

If g n _ x =a(T 1 ,...,T n ;w(v) : v G £(T„_i)), then 



E 



[ J! kj(v)[A {v) + A^X,^ + A 2 (v)X v2 - X v ]I{Vn = v} 
= E[ ^ ro(u)I{F n = u}E[A (?;) + A 1 (u)X„i+A 2 («)X w2 -X„|^ n _i] 

«££(T n _i) 
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and 

E[A (v) + Ai(v)X vl + A 2 {v)X v2 - X v \g n _ x ] = E[Ai + A 2 - 1]E[X ] + E[A ] = 

either in case (b) and (c) or in case (a) when 7 = 1 and E[Xo] = in. Hence, E[VK*] = E[Wg] = 
E[JCo], which completes the proof. □ 

Proof of Proposition \2.S[ Using Proposition 14.11 to the convexity of the Wasserstein distance Z 7 
(7 < 1) and Lemma [5T4l (a) one gets 

^(M*,Moo) < ^^(l-e-'flyCC./ioc) 

n>0 

< e ~ l ( l - e- t )"c„( 7 ); 7 (/2 , Moo) = e-^-^l^Jlo, Moo ). 

□ 

Proof of Proposition \2.4\ Since q(l) < 1, in this case E^i+y^] 7^ 1. Hence, thanks to Proposition 
I4z.ll to (l4"Tj) and to the convexity of the Wasserstein distance 



Hindoo) < E ^i 1 ~ Moo) 

n>0 

< ]T e-*(l - e-*)"c„(l)/i(/io, Moo) = e-' (1 -« (1)) /i(Mo, Moo). 

?!>0 

Furthermore, if E[Xq] = m, then 

I vn t (dv) = Y e ~\ l - e _t ) n E[W*] = Y e _t (l - e"*)"m = m, 

n>0 ri>0 

since E[VK*] = in as stated in (b) of Lemma \h .41 □ 

Proof of Proposition^^ Using Proposition ^. 1[ to the convexity of t G (7 > 1) and Lemma r5T4~l 
(b) one gets 

^(Mt, Moo) < ^ e-*(l - e-*) n q(C, Moo) 

n>0 

< £ e-*(l - e - t )"c„(7)?^(Mo, Moo) = e-^-^Z^Mo, Moo). 



Arguing exactly in as in the previous proof using Lemma 15.41 (b) in place of Lemma 15.41 (a) one 
proves J vfi t {dv) — in. □ 

Proof of Proposition \2~b\ Analogous to the proof of Proposition [2751 using Lemma l5.4l (c) in place 
of Lemma [531(b). □ 
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